Untangling supernova-neutrino oscillations with beta-beam data 
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Recently, we suggested that low-energy beta-beam neutrinos can be very useful for the study of 
supernova neutrino interactions. In this paper, we examine the use of a such experiment for the 
analysis of a supernova neutrino signal. Since supernova neutrinos are oscillating, it is very likely 
that the terrestrial spectrum of supernova neutrinos of a given flavor will not be the same as the 
energy distribution with which these neutrinos were first emitted. We demonstrate the efficacy of the 
proposed method for untangling multiple neutrino spectra. This is an essential feature of any model 
aiming at gaining information about the supernova mechanism, probing proto-neutron star physics, 
and understanding supernova nucleosynthesis, such as the neutrino process and the r-process. We 
also consider the efficacy of different experimental approaches including measurements at multiple 
beam energies and detector configurations. 
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I. INTRODUCTION 

Neutrinos are of fundamental importance during the 
late stages of the evolution of a massive star. Although 
the supernova problem is not yet fully understood, it is 
thought that neutrinos play a crucial role in the core- 
collapse and subsequent explosion. 

When the thermonuclear fuel in the center of the star 
is exhausted, the lack of elements left to burn and pro- 
duce the pressure to maintain its hydrostatic equilibrium 
makes it impossible for the star to prevent the implo- 
sion of the core. Once the core's mass is too large, it 
becomes gravitationally unstable, and starts to collapse. 
The gravitational binding energy that is liberated will 
be released in the form of neutrinos. At first, there is 
a plethora of neutrinos that are generated by the neu- 
tronization processes accompanying the gravitational col- 
lapse. These are free to escape from the star. During 
later stages, the densities and temperatures in the center 
become so high that despite their small interaction cross 
sections, the neutrino diffusion time exceeds the time 
scale of the implosion. The neutrinos are trapped, and 
the equilibrium in the center is extended to weak interac- 
tions. The gravitational binding energy leaks out in the 
form of neutrinos, (which are produced in pairs as well 
as by electron capture), on a neutrino diffusion timescale 
of about ten seconds. A small portion of this energy is 
deposited in the material above the proto-neutron star. 
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New hydrodynamic scenarios [l|, l3| j as well as the inclu- 
sions of energv released by nuclear burning of infalling 
outer shells [3[, can be helpful in creating an explosion. 
Nevertheless, the neutrinos remain an essential energy 
source for the explosion see e.g. [J]. 

Although these neutrinos are only weakly interacting, 
this enormous amount of particles and energy traveling 
through the different layers of the star is able to cause 
a transformation of the elements synthesized during the 
preceding thermonuclear burningjprocesses in the life of 
the star in the neutrino-process [a, |6|. In addition, it is 
thought that a neutrino driven wind occurs at late time, 
which may produce the r-process elements . 

The outcome of the element synthesis in the neutrino 
driven wind, and the prospects for obtaining an r-process 
are quite sensitive to the relative numbers of neutrons 
and protons in the wind. The relative numbers of neu- 
trons and protons in the wind are determined primar- 
ily by the neutrino interactions Vg + n ^>- p + e~ , and 
Ve+p ^ n + e~^ . Small changes in the spectra will trans- 
late into changes in the relative number of neutrons and 
protons and therefore determine the extent to which very 
heavy r-process nuclei can be formed iS, JJ]. Furthermore, 
small changes in the abundance distribution can be af- 
fected by neutrino- nucleus interactions during and at the 
end of the rapid neutron capture process. Because the 
thresholds for these reactions can be quite high, they are 
very sensitive to the neutrino spectra [10, [HI [l3l • The 
same considerations come into play for gamma ray bursts 

[a [mil. 

The neutrino-process in contrast occurs in the outer 
layers of the star, where neutrinos scatter on pre-existing 
nuclei. Certain rare nuclei can be formed by spallation 
processes : a neutrino scatters on a nucleus leaving an 



excited state within the nucleus that then decays by the 
emission of neutrons and/or protons. Therefore it is not 
the capture of neutrinos on free nucleons, but the capture 
of neutrinos on nuclei that is the determine factor for this 
process. Since thresholds are high, changes in the spectra 
due to neutrino oscillations can have a dramatic effect on 
the abundances produced in these processes ^10|. 

As a matter of fact, neutrinos from the next galactic 
core collapse supernova are much anticipated, since they 
are the only particles giving us the chance to see deep into 
the interior of the event and obtain information about 
the processes driving the explosion and the influence of 
neutrinos on the events. Terrestrial supernova-neutrino 
detectors aim at the observation of supernova neutrinos 
through a variety of processes. The need for neutrino- 
nucleus cross sections comes into play. Despite the small 
cross sections, neutrino-nucleus interactions are a pow- 
erful filter for information, owing to their energy, flavor, 
and spin sensitivity [13, 13 ■ 

Currently on line detectors, such as Superkamiokande, 
MiniBoone, and KamLAND are capable of detecting neu- 
trinos from a galactic supernova, for a review see e.g. [l9| . 
For a supernova 10 kiloparsecs away from the earth, a 
heavy water detector, such as the one which operated 
at the Sudbury Neutrino Observatory would see several 
hundred events in all channels, i^g + d —>■ p + p + e~ , 
Ve + d — > n + n + e~ , and ly + d — *■ n + p + v [2C||. 
Superkamiokande will record thousands of events from 
the reaction Ve + p — > n -|- e+ [2l| . In addition there 
will be neutral current events on Oxygen-16, v +^^ O — > 
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detectable by way of the gammas e.g. [2^ |23|. In a pro- 
posed lead detector, such as OMNIS or LAND [Htll], 
the main signals will come from i/e -|-^°* Pb — >208 ^j _|_ g- 
and ly -1-208 PI, ^208 p^ _|_ ^^^ ^^j^j analogue reactions on 
^^Fe. In these reactions the daughter nuclei will typically 
spall one or two nucleons, providing information about 
the energy of the incoming neutrinos. In KamLAND and 
MiniBoone there will be an inverse beta-decay signal on 
protons, as well as a few events from neutrino interactions 
with 12c. 

There are several uncertainties involved in any future 
measurement of supernova neutrinos. These include un- 
certainties in the range of predicted luminosities and 
spectra of the neutrinos, uncertainties in the type and 
degree of neutrino oscillations that occur, and the uncer- 
tainty in the calculation of the detector neutrino-nucleus 
cross sections. Therefore, it is essential to detect neutri- 
nos in all channels : both the neutral current which is 
sensitive to all flavors and the two charged-current chan- 
nels which are sensitive to Ve and De- When the neutri- 
nos are originally emitted, the z/^, Dfj_, v-y and v-y have 
the highest energies, followed by the i/g, and then the 
i/g. Matter-enhanced neutrino transformation will mix 
these spectra differently in the cases of neutrinos and an- 
tineutrinos. Hence, it is important to be sensitive to the 
different neutrino flavors. This provides the only way to 
maximize the information that can be obtained about the 



original spectra. 

The type of neutrino flavor transformation that will oc- 
cur in the supernova is a rapidly developing field. Effort 
is now being directed at understanding fully the neutrino 
background terms [26, 27, 28, 29, 30, 31, 32, 33, 34] and 
their effect on understanding the supernova. Further- 
more, even if these terms are not important, phase effects 
due to multiple resonances come into play [3^, [Sa, l37| . 
This greatly complicates the simple picture of a single 
"H" and "L" resonance [33, [S^- Due to uncertainties 
in neutrino parameters and in supernova conditions, it 
is not yet clear how the neutrinos will transform during 
the relatively early times of the supernova, when most of 
the neutrinos are being released. However, in addition 
to the simplest scenarios of two resonance points in the 
outer layers of the star, realistic possibilities include a 
complete mixing of all flavors [26|, a partial oscillation 
between flavors [23| also called spectral swapping i3^|, 
or a finely grained energy-dependent effect p55!, l3(| . An 
important part of any future observation of supernova 
neutrinos will be to understand this physics. 

The understanding of any astrophysical neutrino- 
nucleus interaction is limited by our understanding of 
its cross section. While the inverse beta decay on pro- 
tons is well understood, there is little data for neutrino- 
nucleus cross sections. The exception is neutrino-nucleus 
cross section data for ^^c [40l . [4ll . WA . l43j measured us- 
ing muon decay at rest (DAR) or decay in flight spectra 
(DIF) spectra where measurements have been reported 
to ^ 10% accuracy. Still, while the exclusive cross sec- 
tions are well understood theoretically, the predictions 
still fail to reproduce some of the inclusive ones (see 
e.g. di, [H, [li, [131). Electron neutrinos from DAR 
are in the same energy regime as supernova neutrinos 
making measurements at the 10% level, such as those 
at the proposed i^-SNS facility [4^ attractive from the 
point of view of calibrating nuclear structure calculations. 
For a few other nuclei, such as Iron, rneasurements have 
been reported at the '^ 40% level [4^. At present, for 
most neutrino-nucleus cross sections, predictions rely on 
shell-model or RPA calculations. In some cases one can 
use other measurements as a guide for the calculation. 
For example charge-exchange data can show the location 
of the Gamow- Teller resonance in the charged-current 
channel, as can the isovector part of the Ml strength in 
the neutral-current channel. Supernova neutrinos how- 
ever are expected to excite also first forbidden transitions 
within the nucleus, for which it is more difhcult to obtain 
useful data from other sources with which to calibrate the 
calculations. 

A recurring issue in the extraction of information from 
neutrino cross section measurements is the absence of 
monochromatic neutrino beams. To some extent, this 
problem can be met by the flexibility offered by beta- 
beam neutrinos. The idea of producing neutrino beams 
from the beta-decay of boosted radioactive ions, stored in 
a storage ring, was first proposed by Zucchelli [50|- The 
main goal of the proposed beta-beam facility was to ex- 



plore the possible existence of CP violation in the lepton 
sector through the comparison of neutrino versus anti- 
neutrino oscillations. This intriguing proposal has trig- 
gered a feasibility study which is now ongoing [5l|, [S^, |53| . 
The sensitivity on the CP violating phase and on the 
third neutrino mixing angle that can be achieved with 
the original as well as other - higher energy - scenar- 
ios is now being explored in great detail (see e.g. [54| . 
[55J and [56J for a review). The energy range interesting 
for core-collapse supernova physics could be covered by 
low energy beta-beams [53| • Several applications of such 
a facility have been studied recently going from nuclear 
structure studie s 157, [58 , 59, 60], the study of fundamen- 
tal interactions [6l|, 62 , l63l . |64 l65l |66| and core-collapse 



supernova physics [57|, |67 1 . 

These applications could exploit either a devoted stor- 
age ring [58| , or low energy neutrino beams at off-axis 



from the standard storage ring [68| . The application dis- 
cussed here particularly benefits from very low ion boosts 
(i.e. 7 = 3,6). Instead of running at such boosts, the 
very low energy neutrino fluxes could be obtained either 
by taking the flux at different parts of a detector at low 
energy beta-beams [60| or with an off-axis detector from 
the standard storage ring |68j . 

In this paper, we examine what can be learned about 
the Galactic supernova neutrinos from the detector sig- 
nal they produce. We follow the approach proposed 
in [67|, where linear combinations of beta-beam spec- 
tra were shown to be able to mimic supernova-neutrino 
energy-distributions very accurately, thus allowing one 
to avoid the problems related to uncertainties and model 
dependencies in theoretical studies, and to the lack of 
monochromatic neutrino beams in experiments. As part 
of our analysis, we explicitly show how it is possible to 
disentangle a supernova-neutrino spectrum which has be- 
come "mixed" due to neutrino oscillations. This can be 
done through a series of carefully chosen beta-beam spec- 
tral measurements. 

The paper is organized as follows : After an introduc- 
tion of the most common descriptions of the supernova 
neutrino spectra and the experimental neutrino energy 
distributions, our technique for the construction of syn- 
thetic spectra is discussed. The quality of the generated 
spectra is evaluated and we suggest some opportunities 
for improvement of the technique in actual implementa- 
tions. Finally, section IIVI elaborates on the information 
about the supernova and about oscillations that can be 
extracted from the neutrino signal in a terrestrial detec- 
tor. 



II. NEUTRINO SPECTRA 

Energy distributions of supernova neutrinos are shaped 
by the circumstances in which the neutrinos are emit- 
ted. Neutrinos leaving the star are responsible for the 
cooling of the proto-neutronstar forming in the star's 
core. Hence, their spectrum resembles a thermal one, 



with temperatures reflecting the conditions at the site 
where they decoupled. However, the fact that different 
kinds of neutrinos are involved in different interactions, 
and that the reactivity of the (anti)neutrino depends on 
its energy, flavor and helicity, modulates this picture. 

For all (anti)neutrino flavors, the energies are in the 
range of a few to a few tens of MeV, although calcula- 
tions of neutrino transport which use different opacities 
achieve somewhat different spectra [70]. Traditionally, 
Fermi-Dirac spectra were put forward as the convenient 
description for the energy distribution nsNi^u^) of su- 
pernova neutrinos 
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Recent results showed the supernova neutrino energy dis- 
tribution to be accurately parametrized with a power-law 
distribution ^,^ 



^iPL[(e),Q](ei^) 



(e.) 



-(a+l).^ 
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where T represents the temperature at the decoupling 
site of the neutrinos and (e^) is the average neutrino 
energy. Up to zeroth order in rj, both quantities are re- 
lated by T = 3.15 (e^). In general, higher average en- 
ergies or temperatures lead to broader spectra reaching 
a maximum at higher energy values, and with enhanced 
tails. The parameters rj and a allow to adjust the width 
w — \/(e^) — (Ei/)^ of the spectrum. Larger values for 
a or rj reduce the width of the spectrum and the influ- 
ence of its tail. Njj is a normalization factor depending 
on the width parameter. The parameterizations allow 
to construct spectra that are equivalent up to their sec- 
ond moment by adjusting the energy and width parame- 
ters. These spectra are illustrated in Fig. [TJ The higher 
energies correspond to mu- and tau-neutrinos, interact- 
ing only through neutral-current interactions and decou- 
pling close to the center of the star at relatively high 
temperatures. The low-energy spectra are important for 
electron neutrinos and antineutrinos, whose opacities are 
built from neutral- and charged-current interactions. The 
precise shape of the spectrum and its tail is very impor- 
tant for the nuclear response to supernova neutrinos, as 
cross sections are rising fast with increasing neutrino en- 
ergies. This is due to the fact that the supernova neu- 
trino energy range is probing the giant resonance region 
of the nuclear spectrum, where cross sections are vary- 
ing fast. This makes nuclei very sensitive probes on one 
hand, but on the other hand makes the concept of nuclei 
as supernova-neutrino detectors very sensitive to uncer- 
tainties in nuclear structure calculations too. 

The experimental neutrino beams mainly used up to 
now to study neutrino interactions at these energies, stem 
from pion decay at rest: 

7r+ ^ /i+ + Vf_,, (3) 



and the subsequent decay of the muon 
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FIG. 1: Neutrino energy-spectra. The upper panels show different power-law parameterizations for supernova neutrinos. Left, 
fixed width and varying average energies (a), right, varying width for two different values of the average energy (b). The lower 
panels show the equivalent Fermi-Dirac distributions. The width w is expressed in units wq — ■^^, for which power-law (with 
a=2) and Fermi-Dirac spectra (77 = 00) coincide. 



The muon neutrinos have energies e^^ — 29.8 MeV, the 
other neutrinos produced in these reactions have energies 
distributed according to the Michel spectrum : 



96£^ 
ml 

324 3 
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providing electron neutrinos and muon antineutrinos 
with energies up to 52.8 MeV. These spectra are illus- 
trated in Fig. [21 Their energy covers the same range 
as that of supernova neutrinos, but the shape is rather 
different. Peak energies are much higher than those for 
supernova neutrinos and the long tail, important for the 
study of supernova neutrinos, is lacking. Fig. [3] shows 
that low-energy beta-beam neutrino energy distributions 
are covering the same energy range as that for supernova 
neutrinos, with a shape that is remarkably similar to that 
of supernova-neutrino spectra. 



III. RECONSTRUCTING NEUTRINO 
SPECTRA AND NUCLEAR RESPONSES 

We aim at optimizing the information that low-energy 
beta-beams can provide about supernova-neutrino inter- 
actions and the interpretation of a supernova-neutrino 
signal in a terrestrial detector To this goal we construct 
normalized linear combinations riM-f of beta-beam spec- 
tra n^. : 
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FIG. 2: Michel spectra, showing the v^ and 77^ energy distri- 
butions stemming from pion decay-at-rest. 
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FIG. 3: Normalized beta-beam neutrino-spectra stemming 
from ^**Ne for different boost factors 7 between 3 and 15, for 
a target with a cross sectional area of 4 m^, located 10 m from 
one end of the ring, which has a straight side length of more 
than 90 m. 



and varied the boost factors 7^=1^... ^tv and the expansion 
coefficients 0^=1^... ^at to minimize the expression 



de^ \ntq^{ev) - nsN{sv)\ 



(10) 



fit 



(e.) 



In this way we obtain a synthetic spectrum 
that is the best fit to the supernova-neutrino energy- 
distribution nsNi^v) for particular values of the aver- 
age energy and width. Using the norm of the expression 
\nN-y{£v) — nsN{£u)\, rather than the square of the norm, 
we avoid giving a larger weight to the peak of the distri- 
bution. This is important in view of the fact that the 
spectrum's tail is more important to our applications. 

Fig. m shows some results of best fits to power-law 
spectra using N=3 and N=5 beta-beam spectra in the 



synthetic energy distribution, for different supernova- 
neutrino energy-distributions. The 7's were allowed to 
assume integer values between 5 and 15. The resulting 
synthetic spectra are compared with the original super- 
nova spectrum, and the corresponding folded cross sec- 
tions 



fold 
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arc confronted. These are especially important because 
the energy dependence of the yield of neutrino interac- 
tions scales directly with the folded cross section. 

The overall fit is quite good. It is clear that, especially 
for spectra with lower average energy, the synthetic spec- 
tra tend to peak at slightly higher energies. This is due to 
the fact that the beta-beam spectrum at 7=5, the lowest 
one included in the fit, is just slightly higher in energy 
than the lowest predictions for supernova neutrino energy 
distributions are. The procedure does not use the free- 
dom to include higher gammas in the constructed spec- 
trum, the fit merely seeks for better agreement at low 
energies. 

It is important to emphasize that the fitting procedure 
is very powerful in reproducing the high-energy tail of 
the supernova neutrino energy spectrum. The inset in 
Fig. m illustrates that, whereas for incoming neutrino en- 
ergies below and around 20 MeV, the synthetic spectra 
are still wobbling around the original one, for energies 
above 30 MeV the fit becomes much more stable. Es- 
pecially for constructed spectra including 5 gamma dis- 
tributions, the agreement is very good. The fact that 
the synthetic spectrum peaks at slightly higher energies 
than the original one does not limit the strength of the 
proposed procedure. The important quantity to repro- 
duce is not the spectrum, but the cross section, folded 
with the spectrum. This folded cross section determines 
the nuclear response to supernova neutrinos. It reaches 
a maximum around 30 to 40 MeV, indicating that neu- 
trinos with these energies have the largest impact on the 
nuclear response [zl] . The bunch of curves right in Fig. 3] 
shows cross sections as a function of the incoming neu- 
trino energy, folded with the supernova neutrino spec- 
tra, and folded with the constructed spectra. The folded 
cross-section curves pop up only in the very tail of the 
energy distribution, and are in very nice agreement. 

Fig. [5] compares total folded cross sections for super- 
nova spectra asN — J de,^ o-{£u)nsN{£v) and the equiva- 
lent synthetic folded cross sections 
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For ^®0, the overall agreement is very good, support- 
ing the strength of the proposed procedure. For the 
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FIG. 4: In each plot, the curves on the left show the original supernova-neutrino spectrum «(<£>, q) (full line), and the fit with 
three beta-beam spectra at different gamma 713^ (dashed), and 5 beta-beam spectra ns-y (dotted). The curves on the right show 
the ^^'Oliyju'Y^O* cross section, folded with the SN neutrino spectrum o-^°>^ (full line), the artificial spectrum constructed 
using 3 different beta-beam spectra 0-3° (dashed), and the artificial spectrum built from 5 beta-beam spectra at different 
gammas CTg° (dotted). The cross section was obtained with a continuum random phase approximation calculation [7J|. The 
inset shows the ratio of constructed spectrum to the original one for the synthetic spectrum containing 3 (dashed line) and 5 
beta-beam spectra (dotted). 



deuteron, the cross sections are much more sensitive to 
lower energies than those for massive nuclei are, and for 
spectra with small average energies the fit overestimates 
the folded cross sections. The agreement becomes better 
between 16 and 18 MeV average energy, growing to an 
almost perfect match for the highest energies we exam- 
ined. 

The fact that the technique is already very efficient in 
the straightforward way it was presented, leaves oppor- 
tunity for optimization. In the following paragraphs, we 
suggest a number of ways to make the technique more 
efficient. 

(1) Figure [S] shows that the folded cross sections be- 
have very smoothly as a function of the boost factor 7. 
This is a major advantage for the proposed procedure : 
a number of measurements at particular gammas will al- 
low to infer other values. Thus, the fitting procedure can 
be offered more freedom in its choice of 7, which will 



definitely result in an even better agreement. Moreover 
the curves in Fig. [S] might be used as guide to extrap- 
olate cross sections to lower gamma values, where the 
event rate in the detector becomes small. All low-energy 
fits, and especially the synthetic spectra for the deuteron, 
would benefit from the inclusion of gammas smaller than 
5 in the fit. 

(2) Once one has experimental beta-beam data avail- 
able, the fitting procedure of Eq. [TO] can be optimized 
by including an energy-dependent weight function in the 
minimization process : 



dsi, \nNy{eu) - nsN'Y{£u)\w{e^). 



(15) 



In this way, it becomes possible to ensure that the ob- 
tained agreement is best in important energy regions i.e. 
the match between original spectrum and fit should be 
closest in energy regions where the cross section is high. 
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FIG. 5: Total cross sections for neutrino scattering on ^'^O (a), and on the deuteron (b), folded with different spectra. The 
^®0 results were taken from Ref. [73|. The folded cross sections are shown for the supernova neutrino spectrum (full line), a 
synthetic spectrum with 3 (dotted line) and with 5 gamma spectra (dashed lines) in the construction. The different bunches 
of curves correspond to different widths of the energy distributions, as indicated by the parameter a. For the deuteron, cross 
sections were taken from Ref. [7g]. The cross sections folded with the original spectra (full line) and with the fit consisting of 
5 gamma spectra (dashed line) are shown for neutral and charged-current neutrino and antineutrino scattering. The width of 
the energy spectrum was kept fixed at q = 3 in this case. 
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FIG. 6: Folded cross section as a function of 7. The left panel shows the results for neutral current neutrino scattering on the 
deuteron j/d— >J^'pn, the right panel for neutral current scattering on Pb. 



As a matter of fact, once experimental values of the 
folded cross section are available, the optimum choice 
for the weight would be a cross section calculation cr*'*, 
translating the minimization procedure in 
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The synthetic spectrum tv^ (si) obtained with this pro- 
cedure serves as a guide for the choice of gammas in 



the experiment. The appropriate combination of the re- 
sponses in the detector then provides a very good repro- 
duction of the supernova-neutrino signal. This is illus- 
trated in Fig. [T] Especially for spectra with smaUer en- 
ergies, using an energy-weighted fit improves the results 
of the procedure. 

(3) Up to now, we have shown that these responses 
can be reproduced using only 3 or 5 /3-beam measure- 
ments at different 7. This restriction is rather aca- 
demic. The smooth behavior of the response as a func- 
tion of the boost factor 7 assures that interpolating ex- 
perimental results would not introduce extra uncertain- 
ties. Hence the proposed procedure will mainly benefit 
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TABLE I: Expansion parameters a^=5_..._i5 for the fit to the power-law supernova-neutrino spectrum, defined by average energy 
(e) and width a. The quantity J de,/|n-'*'(e^) — nsN{£v)\ in the last column provides a measure for the goodness of fit. 
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TABLE IL Expansion parameters for the fit to supernova-neutrino spectra, using the fluxes from Ref. [69|, obtained at different 
parts of a detector placed 10 m away from a 450 m storage ring with a straight section of 150 m. All considered parts of the 
detector are cylindrical and have a volume of 71.2 m'^, radius and thickness of the parts are given beteween brackets. The 
quantity / de^|n^'*(£i/) — nsN{£v)\ in the last column provides a measure for the goodness of flt. 



from the advantages using a larger set of expansion func- 
tions brings along. Table U shows some results of best 
fits to power-law spectra, for different supernova-neutrino 
energy-distributions, where all 7 values between 5 and 15 
were allowed to contribute to the synthetic spectra. As 
mentioned above, the most important quantity to study 
is the differential folded cross section i.e. the folded cross 
section as a function of the final energy of the target 
nucleus. This provides a measure for the energy trans- 
fer from the neutrinos to the material they are inter- 
acting with, and it determines the reaction products in 
nucleosynthesis processes. For supernova neutrino detec- 
tion, the differential folded cross sections indicate what 
one will 'see' in the detector, as the energy transfer and 
the excitation energy of the nucleus determine the decay 
products that will be observed in the detector. Fig. [5] 
illustrates these response in terms of differential cross 
sections for neutral-current scattering off ^"^Pb, show- 
ing an almost perfect agreement between original and 
fitted response where all 7- values between 5 and 15 con- 
tribute to the synthetic spectra. As a matter of fact, the 
figure predicts the supernova-neutrino signal in a terres- 
trial detector and the response on nuclei of interest in the 
neutrino nucleosynthesis processes in terms of beta-beam 
responses. 

(4) Table U shows that the fit to the lowest-energy 
supernova neutrino spectra is dominated by beta-beam 



spectra at very low gammas. From an experimental point 
of view, it might be beneficial to use fluxes obtained at 
an off-axis detector or at different parts of the detector as 
an effective way to lower the average energy of the neu- 
trinos, rather than running a facility at very low boost 
factors. Table [TTl uses the fluxes from Ref. [6^ obtained 
at different parts of a detector to vary the average energy 
and shape of the synthetic spectra, and shows that this 
is an excellent way to obtain a fit with the same accuracy 
level for a beta-beam facility run at higher gammas. 



IV. TERRESTRIAL DETECTION OF 

SUPERNOVA NEUTRINOS AND NEUTRINO 

OSCILLATIONS 

The very satisfying agreement between original and 
constructed spectrum, especially for massive nuclei, sug- 
gests that it is possible to construct supernova-neutrino 
signals and neutrinonuclcosynthesis cross sections, sim- 
ply by taking the results of the beta-beam measurement 
in the detector, without going through the intermedi- 
ate step of using the nuclear structure calculation. The 
procedure is easy. For each set of beta-beam data at a 
given gamma, there will be a measured response in the 
detector. For neutral current on ^^O and ^°*Pb it will be 
spalled neutrons, protons and coincident gammas, for the 
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FIG. 7: Folded cumulative cross section for neutral current 
neutrino scattering off ^^0 as a function of the excitation 
energy uj. The figure compares cross sections folded with dif- 
ferent spectra : original power law spectra (full line) , the syn- 
thetic spectrum obtained without weight in the fit (dotted), 
and the fit with the cross section calculation used as weight 
(dashed line). 



charged-current channels it will be electrons or positrons 
plus possibly coincident neutrons, protons and electrons. 
Our results for the differential cross sections now clearly 
show that taking the appropriate linear combinations of 
the measured response provides a very accurate picture of 
the neutrino-nucleosynthesis processes or of the response 
of the detector to an incoming supernova-neutrino spec- 
trum. 

The accuracy reached by the fit ensures that it be- 
comes possible to discriminate between supernova neu- 
trino spectra that are quite similar [63|- Fig. [§] illustrates 
that even for equivalent power-law and Fermi-Dirac spec- 
tra, identical up to the second moment of the distribu- 
tion, the fit allows to distinguish between both responses. 
The procedure is sensitive to very subtle differences in the 
spectrum. 

In fact, as was shown in Ref. [67| the discriminative 
power of the procedure is such that it becomes possi- 
ble to reconstruct the incoming supernova neutrino spec- 
trum from the response in a terrestrial neutrino detector. 
Hence, this technique provides a way to study the super- 
novaneutrinos without having to rely on the intermediate 
step of theoretical cross section estimates with their re- 
lated uncertainties. 

This opens ways to a model-independent study of os- 
cillations in supernovaneutrinos. In terms of energy dis- 
tributions, the supernova neutrinos are emitted with a 
low-energy and a high-energy component, the first asso- 
ciated with electron (anti)neutrinos, the latter with the 
heavy-flavor mu- and tau neutrinos. In a terrestrial de- 
tector, both types will interact through neutral-current 
reactions. Charged-current detection channels are only 



accessible for electron-type neutrinos, as the production 
of a massive mu- or taulepton requests more energy than 
is available for supernova neutrinos. If we denote the low- 
energy components by n''a^{e,y^), n^sNi^v^) and the high- 
energy components by n|^ (^i/,, , £i. J , and n^^ [e^^ , ei7^ ) , 
the original energy distributions can be written as 

nsN{£u) = n^sN{£^.)+'^n%j^{e^^^,e^^), (18) 
nsN{£-) = nl;^(£i7j-H 2n^^(ei7^,£T7j. (19) 

If we assume that neutrinos of all flavors are emitted in 
equal numbers, there will be twice as much high-energy 
neutrinos. Oscillations will transform the high and low- 
energy components in a different way. This does not 
affect the neutral current signal, as it is not sensitive to 
the neutrino flavor. In terms of the energy spectra, the 
neutral current signal can then be written as 

Snc{^)^ dsy (nsjv(e,.J -t- 2n|^(£^^,£^J 

x <TNc{£^,^), (20) 

where the energy integration runs over all neutrino fla- 
vors. The energy distribution of the neutrinos responsible 
for the charged current signal will be distorted according 
to : 

(21) 
with R denoting the fraction of electron-type neutrinos 
that remained unchanged and / the fraction of heavy- 
flavor neutrinos that were transformed into electron neu- 
trinos. In the equation above R + I = 1, for any given 
energy of the neutrinos. For simplicity, in our example 
below, we assume that the neutrinos will evolve either 
completely adiabatically or completely non-adiabatically 
so that in addition, both R and / are constant as a func- 
tion of neutrino energy, although our analysis could be 
expanded to include the more general case. The param- 
eters R and / contain the information about oscillations 
that the supernovaneutrinos are carrying and can be re- 
covered from the signal in the detector in two steps. First, 
from a set of constructed spectra unj, the two linear 
combinations of beta-beam responses 
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FIG. 8: Comparison between differential cross sections for neutral current scattering on Pb, folded with a power-law 
supernova-neutrino spectrum (full line) and synthetic spectra taken from Tab. |T] (dashed line) with /3-beam components between 
7=5 and 15, for different energy distributions : (£i=14 MeV, a=3 (a), (e)=22 MeV, a=3 (b), {e)=18 MeV, q=2 (c), and {£)=18 
MeV, Q=4 (d). Cross sections were taken from [74| ]. 



minimizes the difference with the neutral current detec- 
tor signal from Eq.l^Dlare selected. In the neutral-current 
signal, we have not distinguished the antineutrino compo- 
nent from the neutrino component. The antineutrino sig- 
nal is not too different from its neutrino counterpart, and 
so we have disregarded these differences in this illustra- 
tion of the proposed reconstruction method. The result- 



ing two sets of expansion parameters (a: 



fit, I . ,■ _ 



and (af 
nations 



fit,h ■ 



3,15) 



3, 15) are then used to construct combi- 



de,(i?^^* n'^/"(e,) + I-^^* n'*'^^*(e.))fTcc(e.), (25) 



and seek for the oscillation parameters i?-''** and /■''** that 
yield the best agreement between the expression of Ea.[25l 
and the signal of Eq. [51] 

The results of this procedure are shown in Fig. 1101 
where the parameters R and / were extracted from a 
supernova- neutrino signal consisting of two mixed power- 
law distributions. The reconstruction was realized using 



a library of 125 spectra nJ^Ls 15 (^i^) with average en- 
ergies between 12 and 24 MeV and width varying from 
a = 2 to 5. The method is very successful, even when 
the spectra are very close together as illustrated in panel 
(c) of the figure. 

As for most neutrino experiments, especially at rela- 
tively low energies, statistics can be expected to be rather 
poor. This will result in uncertainties on the observed 
supernova-neutrino signal and on the beta-beam data. 
We investigated the influence of these uncertainties on 
the reconstruction procedure by adding random noise 
to the 'data points' of the signal of Eqs. [201 and [51] i.e 
Snc{^) and Scci'-^) were replaced by S{uj) * {l.+p*r), 
where r was chosen radomly between and 1 for each 
data point and p determines the size of the noise, p = 0.05 
and p = 0.10 for the simulations shown in Fig. [11] . The 
accuracy of the first step in the procedure , i.e. the re- 
construction of the spectra n^sNi^:^) and ngj^{e,^) from 
the neutral current signal Snc{(^) is of course crucial for 
the subsequent determination of R and /. It was already 
shown in Ref. [63| that the procedure for reconstruct- 
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FIG. 9: Cumulative differential folded cross section for scat- 
tering off Pb. The figure compares power-law (full line) 
and Fermi-Dirac spectra (dashed) with the same average en- 
ergies and widths and their respective fits (dotted for the fit to 
the power-law spectrum, dash-dotted for the fit to the Fermi- 
Dirac spectra). The parameter w refers to the width relative 
to that of a Maxwell-Boltzmann distribution with 10=1. 



Concluding, we have shown that beta-beam measure- 
ments can be of great help in the analysis of a supernova- 
neutrino signal. Making use of the information provided 
by beta-beam data provides a model-independent way for 
the extraction of information on the supernova-neutrino 
energy-distribution from its signal in a terrestrial detec- 
tor. As neutrinos are the only particles reaching us from 
the core of the exploding star, this makes beta-beams 
an important instrument for probing the processes driv- 
ing the core-collapse supernova and the circumstances 
reigning in the heart of the event. The central idea of 
the proposed technique consists of a reconstruction of 
supernova-neutrino responses using experimental beta- 
beam data, hence it is also of direct use for the prediction 
of neutrinonucleosythesis reaction products. We propose 
to use this technique as a lever to make advantage of beta- 
beam data in disentangling the oscillation characteristics 
from mixed spectra. We have shown that the procedure 
is effective in extracting oscillation parameters from the 
signal in a terrestrial detector, illustrating the use of the 
technique for gathering information about the supernova 
as well as its neutrinos. 



ing a single spectrum is stable against uncertainties at 
this level. Fig. [Tl] corroborates these results and shows 
that the second step, the reconstruction of the parame- 
ters R and / from the fit to the charged current spec- 
trum is also stable against these uncertainties. For a fit 
to 100 different spectra, all obtained by adding random 
noise to the signal, the results for i?^** are strongly con- 
centrated around the original i?-values. In agreement 
with the results of Ref . [63| , the reconstructive power of 
the proposed procedure slightly dwindles for spectra at 
lower energies where the influence of uncertainties be- 
comes larger. 
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FIG. 10: Extraction of the oscillation parameters J?'*' and /■''* from a fit of the expressions in Eqs. (|24|) and (|25|) to 
the neutral Snc{(^) = J^ de„{nsff{£i,) + 2nspf{£„))aNc{£v,<^), and charged-current Sccii^) = J^ d£„{R n'sMi^i^) + 
I nsN{sv))'^cc{£v,(j-') supernova-neutrino signal in an oxygen detector. For each fitted point the original parameter values 
{R, I) are given between brackets. The agreement between {R, I) and the reconstructed (i?-'**, J-*'*) is very good. 
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FIG. 11: Influence of noise on the reconstruction of the pa- 
rameter R from the supernova neutrino signal in a terrestrial 
detector. The figure shows the value of W'^, for original 
values _R=0.2, 0.5 and 0.7 for 100 fits to a signal subject to 
random noise at 5 and 10% level, as described in the text. 
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